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The existence of periodic orbit bunches is proven for the diamagnetic Kepler problem. Members 
of each bunch are reconnected differently at self-encounters in phase space but have nearly equal 
classical action and stability parameters. Orbits can be grouped already on the level of the symbolic 
dynamics by application of appropriate reconnection rules to the symbolic code in the ternary alpha- 
bet. The periodic orbit bunches can significantly improve the efficiency of semiclassical quantization 
methods for classically chaotic systems, which suffer from the exponential proliferation of orbits. 
For the diamagnetic hydrogen atom the use of one or few representatives of a periodic orbit bunch 
in Gutzwiller's trace formula allows for the computation of semiclassical spectra with a classical 
data set reduced by up to a factor of 20. 
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I. INTRODUCTION 

Semiclassical theories provide the link between quan- 
tum spectra and the dynamics of the underlying classi- 
cal system. For systems with chaotic classical dynam- 
ics Gutzwiller's trace formula [1] expresses the density 
of states in terms of parameters of the isolated unsta- 
ble periodic orbits. Chaotic systems typically exhibit 
an exponential proliferation of the number of periodic 
orbits with growing length. The periodic-orbit sum of 
Gutzwiller's trace formula usually does not converge, and 
special techniques as, for example, the cycle expansion 
[2] or harmonic inversion [3]-[5] method must be applied 
for the computation of semiclassical eigenvalues from a fi- 
nite set of trajectories. Although these techniques can be 
very powerful for certain systems a remaining conceptu- 
ally weak point of periodic orbit quantization is still that, 
due to the exponential proliferation of orbits, the gradual 
improvement or extension of semiclassical spectra usually 
requires an exponentially increasing set of periodic orbit 
data. In this article we propose a way to significantly 
reduce the number of trajectories which is necessary for 
periodic orbit quantization by employing properties of 
the orbits which have recently been investigated in con- 
nection with universal features of quantum chaos. 

According to the Bohigas-Giannoni-Schmit conjecture 
[6] a signature of quantum chaos is universal spectral 
fluctuations on the scale of the mean level spacing. For 
example, for systems with time reversal symmetry the 
spacings show a Wigner distribution. On the way to- 
wards an understanding or even proof of the conjecture 
the correlations between orbits which play an important 
role in the semiclassical approximation to the spectral 
form factor have been studied [7HT5]. Strong correla- 
tions exist only between orbits which have identical or 
near identical actions. It has been found that periodic 
orbits consist of different segments [14] . In each seg- 
ment, an orbit follows closely its neighboring orbit or 
the time-reverse of this orbit, but the orbits differ in how 
these segments are connected at self-encounters. As a 



consequence long periodic orbits of hyperbolic systems 
do not exist as independent individuals but rather come 
in closely packed near action-degenerate periodic-orbit 
bunches. Pairs of orbits with two differently connected 
loops provide the leading off-diagonal contributions to 
the form factor [8]. The role of self-encounters and 
periodic-orbit bunches for universal level correlations in 
quantum chaos has been investigated in detail in Refs. 

EM. 

Here we want to demonstrate that near action- 
degenerate periodic-orbit bunches can help to signifi- 
cantly improve the efficiency of semiclassical quantiza- 
tion methods. The idea is not to use all individual pe- 
riodic orbits up to a given length but only one or very 
few representatives of a periodic-orbit bunch. Those rep- 
resentatives are appropriately weighted in the periodic 
orbit sum according to the size of the bunch, which can 
be determined, for systems with a known symbolic code, 
solely on the level of the symbolic dynamics without a 
numerically expensive periodic-orbit search. 

Results are presented for the hydrogen atom in a mag- 
netic field, which has a long history as a prototype model 
of a real physical system with signatures of quantum 
chaos [T5J [TO] . At sufficiently high positive energies it is 
an open system with fully hyperbolic classical dynamics, 
and a unique symbolic code for the periodic orbits does 
exist [P7HT5] . Semiclassical resonances of the diamag- 
netic hydrogen atom have already been obtained with a 
modified and extended cycle expansion technique 20 . 

The article is organized as follows. In Sec.|n]we discuss 
the classical dynamics of the diamagnetic hydrogen atom 
and present an efficient multi-shooting algorithm for the 
numerical periodic-orbit search. In Sec. |III| we reveal the 
existence of periodic-orbit bunches and introduce four 
reconnection rules which allow for the grouping of or- 
bits on the level of the symbolic dynamics. Semiclassical 
resonances are computed in Sec. |IV| with the harmonic 
inversion method, and it is illustrated how using periodic- 
orbit bunches significantly improves the efficiency of the 
method. Concluding remarks are given in Sec. [Vj 



2 



II. CLASSICAL DYNAMICS AND 
PERIODIC-ORBIT SEARCH 

In this section we recapitulate the basic equations for 
the classical dynamics of the diamagnetic hydrogen atom 
and present a very efficient multi-shooting algorithm for 
the periodic orbit search. 

In atomic units the Hamiltonian of the diamagnetic 
hydrogen atom reads 

H= 1 -p*-l + 1 - 1 L z + 1 - 1 *(x 2 + y 2 )=E > (1) 

where 7 = S/(2.35 x 10 5 T) with B the magnetic field 
strength. We only consider states with magnetic quan- 
tum number m = L z = 0. The classical dynamics 
does not depend separately on the energy E and the 
magnetic field strength 7 but only on the scaled energy 
E = E7- 2 / 3 . Using a regularization of the Coulomb sin- 
gularity in scaled semiparabolical coordinates [El US] 

H = ~/ 1/3 \/r + z and v = ~/ 1/3 y/r - z , (2) 

the classical equations of motion in the transformed time 
r, with dr = 2rdt, are derived from the scaled and regu- 
larized Hamiltonian 

h=\{pl+pl)+V{pL,v) = 2, (3) 

with the potential 

V(jm, v) = —E + v 2 ) + \ (A 2 + A 4 ) • (4) 

The numerical computation of periodic orbits is a pre- 
requisite for the semiclassical quantization of chaotic sys- 
tems. In chaotic billiards with an existing symbolic dy- 
namics such as the three-disk scatterer or the hyperbola 
billiard, a periodic orbit with given symbolic code can be 
computed efficiently by moving the reflection points until 
the orbit length becomes a minimum. For systems with 
smooth potentials the periodic orbit search is more diffi- 
cult. Using a simple shooting algorithm trajectories with 
varying starting points are integrated numerically until 
the initial and final point match. When applying this 
method for a systematic periodic orbit search some or- 
bits are typically found many times while others may be 
overlooked. In particular, long and very unstable orbits 
are hard to find. 

For the computation of the periodic orbits of the dia- 
magnetic hydrogen atom we employ a multi-shooting al- 
gorithm which is adapted to the symbolic dynamics of 
the system. With this algorithm we can find selectively 
each orbit corresponding to a given symbolic code, even 
when the orbit is very long and unstable. 

If one allows the semiparabolic coordinates (/x, v) to be 
positive or negative, the potential |4]) has a C^ v symme- 
try. Periodic orbits can be described by a ternary sym- 
bolic code |17) in a similar way as for the hyperbola or 
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FIG. 1. (Color online) In the multi-shooting algorithm each 
orbit is split into short segments starting and ending on one 
of the axes. The concatenation of various types of segments 
is guided by the symbolic code of the orbit. 



four-disk billiard with the same symmetry. At high ener- 
gies E > E c — 0.329 there is a one-to-one correspondence 
between orbits and the symbolic code, whereas below the 
critical energy orbits undergo bifurcations and the sym- 
bolic dynamics is pruned [T5J [T15] . 

The multi-shooting algorithm for the periodic orbit 
search basically works as follows. Each periodic orbit is 
split into segments, with the number of segments equal 
to the string length L of the symbolic code. The seg- 
ments start on either the /x or the v axis and end on one 
of the axes as illustrated in Fig. [T] When we start all seg- 
ments with initial guesses chosen in accordance with the 
symbolic code the orbit is discontinuous between the end 
point of one and the starting point of the next segment. 
Now the initial conditions of all segments are iteratively 
changed to remove the discontinuities. This problem can 
be formulated as a 2L-dimensional root search. (Note 
that two parameters can be changed for each segment to 
vary the initial conditions on one of the axes.) 

Because the segments are short the multi-shooting al- 
gorithm usually converges rapidly and uniquely to the 
periodic orbit selected by the symbolic code, provided 
the orbit is not pruned or very close to bifurcation. The 
method works very well even for very long and unstable 
orbits. 

An example of a periodic orbit at scaled energy E = 0.5 
is presented in Fig. [2j For the symbolic code of a pe- 
riodic orbit we adopt the ternary alphabet of Ref. [T7] 
with a slightly modified way of the symmetry reduc- 
tion. In the full coordinate space [see Fig. [2]Ja)] where 
the semiparabolic coordinates are extended to negative 
values each orbit is labeled by a sequence of numbers 1 to 
4 for the reflections at the potential barrier in one of the 
quadrants. The C^ v symmetry of the potential allows for 
symmetry reduction of the symbolic code. Symbols zx, 

and are used for the orbit crossing the ix- or ^-axis 
in a clockwise or anticlockwise turn, or moving to the 
opposite potential barrier with a crossing of both axis. 
At this level of symmetry reduction orbits are located in 
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FIG. 2. (Color online) Example of a periodic orbit displayed in (a) the full coordinate plane, (b) the 1st quadrant, and (c) the 
fundamental domain. 



the first quadrant [see Fig.[2^b)]. Finally, an axis symbol 
fi or v is replaced with a + or - sign when it follows a 
different or the same axis symbol (ignoring any interme- 
diate symbols), respectively. This alphabet describes 
the symmetry reduced orbit in the fundamental domain 
[see Fig. [2jc)]. The various levels of symmetry reduction 
are helpful to illustrate the periodic orbit bunches and 
reconnection rules in Sec. IIII1 



III. PERIODIC-ORBIT BUNCHES 

In this section we first provide an example of a group 
of near action-degenerate orbits of the diamagnetic hy- 
drogen atom. We then present four reconnection rules 
which are based on the symbolic dynamics of the peri- 
odic orbits and can be used for grouping the trajectories 
in periodic-orbit bunches. The properties of the periodic- 
orbit bunches are discussed. 




FIG. 3. (Color online) Example of a periodic-orbit bunch 
consisting of 16 near action-degenerate orbits drawn in the 
fundamental domain of the (fi, coordinate space for scaled 
energy E = 0.5. The similarities between orbits are obvious 
in the overview. 



A. Example of a near action-degenerate 
periodic-orbit bunch 

We start our discussion of periodic-orbit bunches of the 
diamagnetic hydrogen atom by way of example of a group 
of 16 orbits with cycle length L = 10 at scaled energy 
E = 0.5. The orbits are presented in Fig. [3] in the funda- 
mental domain of the coordinate space. In the overview 
the similarities between the orbits are quite obvious. All 
orbits appear to run nearly parallel and thus to be located 
in the same area of the phase space. For a more detailed 
discussion the 16 orbits of the periodic-orbit bunch are 
drawn separately in Fig.|4]in semiparabolic (fi, v) coordi- 
nates without symmetry reduction along the fi = v sym- 
metry line. Each orbit is labeled by the symbolic code 
using the ternary alphabet introduced in Sec. [TlJ To make 
the graphs concise only one fundamental period is shown 
in Fig. [4j i.e., some orbits are periodic only in the fully 
symmetry reduced fundamental domain. At first glance 



the orbits in Fig.[4]do not appear to be very similar, how- 
ever, the similarities are obvious in Fig. [3] after complete 
symmetry reduction in the fundamental domain of the 
coordinate space. A more detailed inspection of Fig. [4] 
allows us to reveal both the similarities and differences of 
the members of the periodic-orbit bunch and to illustrate 
the reconnection rules. 

The basic mechanism for the formation of near action- 
degenerate periodic-orbit bunches is that each periodic 
orbit can be split into segments where the starting and 
end points of segments can be connected in various ways 
using simple combinatorial rules. Two types of segments 
exist where the trajectory is located either close to or 
far away from the /i = v symmetry line. They can be 
described most efficiently by the symbolic code of the or- 
bits. A sequence of one or more zero symbols followed by 
one plus or minus symbol characterizes a segment located 
near the fi = v symmetry line, a sequence of plus and mi- 
nus symbols (without any zero) characterizes a segment 



4 



(a) 00--0--+-- (b)00--0+-+- ,. (c)OO--O— +- (d) 00— +-0-- 




(e)00— +-0+- , ' ' (f)00+-0— +- ,. (g) 00+-0--+-- (h) 00--+-0-- 




(i) 00+-+--0+- G)00+-0+-+- . (k) 00+-+--0-- (I) 00-+--0+- 




(m) 00--0+-+- , (n) 00+-0+--+- ,. (o) 00+-+-0-- ,.■*' (p) 00+-+-0+- 




FIG. 4. (Color online) Example of a near action-degenerate periodic-orbit bunch at scaled energy E = 0.5 consisting of 16 
orbits with cycle length L = 10. The trajectories are drawn in semiparabolic (fi, v) coordinates and are labeled by the symbolic 
code. See text for discussion. 



away from the symmetry line. All orbits in Fig. [4 thus 
consist of four segments, e.g., 00-|-|0-|-+ — in Fig. 4ta). 
The orbits in Fig. |4|b)-(p) are obtained from the orbit 
in Fig. |4ja) by reconnecting the four segments in var- 
ious ways. For example, in the orbit 00-|-+ — 1 — | — in 
Fig. |4jh) the two segments - and -+ — are interchanged, 
and in the orbit 00-|-|0-| — +- in Fig. |4^c) the segment 
-+ — is traversed backward. Furthermore, the plus or 
minus symbol following a sequence of zero symbols can 
be interchanged [see orbits 00-|-|0+|-+ — , 00+|-|0-|-+ — , 



and 00+|-|0+|-+— in Fig.gb), (g), and (j)] which is re- 
lated to the mirroring of segments at the [i — v symme- 
try line. The remaining orbits in Fig. [4] are obtained by 
combined application of several reconnection rules, e.g., 
in Fig. [4jk) the segments - and -+ — are interchanged 
and the first segment 00- is replaced with 00+. 
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B. Reconnection rules 

We denote a sequence of one or more zero symbols fol- 
lowed by one plus or minus symbol as (O)-stretch, and 
an arbitrary sequence of plus and minus symbols without 
any zero symbol as (+-)-stretch. The grouping of trajec- 
tories of the diamagnetic hydrogen atom in periodic-orbit 
bunches can now be described by the following four re- 
connection rules. 

Rule R-1: The last symbol of a (0) -stretch, which is 
a plus or minus symbol, can be replaced with the other 
one. 

Rule R-2: The order of symbols in a (+-)-stretch can 
be reversed. 

Rule R-3: In a symbolic code with two or more (0)- 
stretches different (0) -stretches can be interchanged. 

Rule R~4: In a symbolic code with two or more (+-)- 
stretches different (+-) -stretches can be interchanged. 

To illustrate the reconnection rules we discuss once 
more the periodic orbits shown in Fig. [4] The first orbit 
00 — — + — in Fig.Qa) can be decomposed into the (0)- 
stretch 00-, the (+-)-stretch -, the (O)-stretch 0-, and 
the (+-)-stretch -+ — . The remaining 15 orbits in Fig. [4] 
can be obtained from that orbit by application of one or 
more of the reconnection rules. For example, the orbit in 
Fig. [4jb) is obtained by applying rule R-1 to the second 
(0) -stretch 0-, the orbit in Fig. [4^c) is obtained by ap- 
plying rule R-2 to the (+-) -stretch -+ — , and the orbit in 
Fig. |4|d) is obtained by applying the rules R-2 and R-3, 
i.e., the (+-)-stretch -+ — is reversed and then the two 
(+-)-stretches — +- and - are interchanged. Applica- 
tions of the reconnection rules mean that parts of the or- 
bit drawn in semiparabolic (p, v) coordinates are approx- 
imately mirrored at the /i = v symmetry line, parts of 
the orbit are traversed backward, or parts of the orbit are 
visited in a different order. However, in the fundamental 
domain of the (fi, v) coordinate space the global shape of 
orbits belonging to the same periodic-orbit bunch is very 
similar, as can be clearly seen in Fig. [3] 

The reconnection rules are schematically illustrated in 
Fig. [5] A periodic-orbit bunch is characterized by a set of 
k (+-)-stretches, which is a sequence of plus and minus 
symbols (an empty set is allowed), and (for k > 2) the 
same number of (O)-stretches. The (O)-stretches are lo- 
cated in the vicinity of the /1 = v symmetry line and are 
indicated by red solid lines in Fig. [5] The (+-) -stretches 
are located in regions with /i^y away from the symme- 
try line and are indicated by green dashed lines. Each 
(O)-stretch is followed by a (+-)-stretch and vice versa, 
the connections are indicated by the thin black lines in 
Fig. [5] The individual members of the periodic-orbit 
bunch differ in the way how the stretches are connected: 
(+-)-stretches can be mirrored at the /i = v symmetry 
line (rule R-1), traversed backward (rule R-2), and the 
orders of the (0)- or (+-) -stretches can be interchanged 
(rules R-3 and R-4). 

How many orbits belong to a certain periodic-orbit 
bunch? Although we cannot give a precise answer it is 
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FIG. 5. (Color online) Sketch of the four reconnection rules. 
In the coordinate space (O)-stretches indicated by red 

solid lines are located in the vicinity of the fi = v symmetry 
line, and (+-) -stretches indicated by green dashed lines are 
located in regions with (1 / !/ away from the symmetry line. 
The connections between (0)- and (+-) -stretches are drawn 
with thin black lines. (+-) -stretches can be mirrored at the 
(X = v symmetry line (rule R-1), traversed backward (rule 
R-2), and the orders of the (0)- or (+-) -stretches can be 
interchanged (rules R-3 and R-4). 

easy to estimate an upper bound of that number from 
the combinatorial reconnection rules. Let again k be the 
number of (+-)- and (O)-stretches in a symbolic code. 
According to rule R-1 the last symbol of a (O)-stretch, 
which is + or -, can be replaced with the other one, re- 
sulting in 2 k possibilities. According to rule R-2 the or- 
der of symbols in a (+-) -stretch can be reversed. Each 
(+-)-stretch which is not symmetric under that trans- 
form provides a factor of two, resulting in at most 2 fc 
possibilities. The interchange of (O)-stretches and (+-)- 
stretches according to rules R-3 and R-4 provide at most 
(substrings may be identical) kl and (k— 1)! possibilities, 
respectively, where, in the second factorial, it has been 
taken into account that a cyclic permutation of all (0)- 
and (+-) -stretches does not change the periodic orbit. 
Altogether we obtain a maximum number of orbits 

7V™ ax = fc!(/c-l)!2 2fe (5) 

building one periodic-orbit bunch. The 16 orbits shown 
in Fig. |4]have k = 2 (O)-stretches. The number of orbits 
is less than the maximum iV™ ax = 32 because the (+-)- 
stretch consisting of a single minus symbol coincides with 
its reverse. The maximum number iV™ ax increases very 
rapidly with k, e.g., for k — 4 bunches with up to iV™ ax = 
36864 orbits exist. However, for such a huge bunch the 
four (0) -stretches must have different lengths and the 
four (+-) -stretches including their time reversals must 
be different, which is found to be possible only for very 
long orbits with code lengths L > 26. 
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FIG. 6. (Color online) Distribution histogram of the classical 
actions for the 16 trajectories of the near action-degenerate 
periodic-orbit bunch shown in Figs. [3] and [4] The small differ- 
ences between the actions of individual orbits are only visible 
in the enlarged inset. 



C. Properties of periodic-orbit bunches 

The most important feature of the periodic-orbit 
bunches is that the individual trajectories have similar 
shape in coordinate space (see Fig. |3|, and thus all or- 
bits have nearly the same action and stability properties. 
For the 16 orbits shown in Figs. [3] and [4] the near de- 
generacy of the actions is illustrated in Fig. [6j In low 
resolution the distribution of the actions exhibits a sin- 
gle peak at s/2ir « 11.1. Only in high resolution (see the 
inset in Fig. [6]) the distribution of the actions in a small 
but nonzero range As can be observed. 

For a given periodic-orbit bunch the action range As 
basically depends on the number of zero symbols in the 
shortest (0) -stretch. We illustrate this for the orbit pairs 



0000- — — — and 0000+; 

n L—n—1 n L — n— 1 

which consist of a (0) -stretch with n zero symbols fol- 
lowed by a plus or minus symbol and a (+-) -stretch with 
L — n—1 minus symbols. The total length of the sym- 
bolic code is L. For fixed L and n the two orbits belong to 
the same periodic-orbit bunch (see rule R-l). The action 
difference As between the two orbits for various values 
of L and n is shown in Fig. [7] Evidently, As decreases 
exponentially with the number n of consecutive zero sym- 
bols but only weakly depends on the total code length L. 
This result does not depend on the special choice of the 
(+-)-stretch as a sequence of minus symbols, i.e., simi- 
lar results have been obtained for other types of (+-)- 
stretches. A more detailed investigation shows that the 
similarities of orbits in a periodic-orbit bunch become 
more and more pronounced with increasing length of the 
shortest (O)-stretch. The reason is the decrease of the 
angle or distance of the orbit at self-encounters with in- 
creasing length of the (O)-stretches [2"T] . 




FIG. 7. (Color online) Action difference As between two or- 
bits consisting of a (0) -stretch with n zero symbols and a 
(+-)-stretch with L — n — 1 minus symbols. Evidently, As 
decreases exponentially with the length of the (0) -stretch. 



IV. SEMICLASSICAL QUANTIZATION WITH 
PERIODIC-ORBIT BUNCHES 



As already mentioned in Sec. [T] self-encounters and 
periodic-orbit bunches play an important role for the un- 
derstanding of the level statistics of quantum systems 
with an underlying chaotic classical dynamics. We will 
now demonstrate the relevance of periodic-orbit bunches 
for the semiclassical quantization of the diamagnetic hy- 
drogen atom. 

The energy eigenvalues of integrable classical systems 
can be obtained by semiclassical torus quantization. In 
chaotic systems torus structures in phase space are ab- 
sent. It has been shown by Gutzwiller |T] that the unsta- 
ble periodic orbits are the skeleton for the semiclassical 
quantization of chaotic systems. However, the problems 
with Gutzwiller's trace formula are twofold: First, the 
periodic orbit sum is not absolutely convergent and spe- 
cial techniques such as cycle expansion [2] or harmonic 
inversion [3]-[5] methods are necessary to obtain converged 
results. Second, in chaotic systems the number of peri- 
odic orbits grows exponentially with the orbital lengths. 
It is very unsatisfactory that an exponentially increasing 
set of classical data information is necessary to resolve a 
few more semiclassical eigenvalues. However, as already 
mentioned, the required classical data can be significantly 
reduced when using the periodic-orbit bunches instead of 
the individual periodic orbits for the semiclassical quan- 
tization. 



A. Harmonic inversion method 

We employ the harmonic inversion method for periodic 
orbit quantization [3-5, 22J, and, for the convenience of 
the reader, here briefly recapitulate its basic ingredients. 
According to [T] the semiclassical response function of 
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chaotic systems is given by 



po 



A po (E) e lb - 



(6) 



where g® c is a smooth function of the energy, S po is the 
classical action of a periodic orbit, and A po is the am- 
plitude of that orbit including phase information. The 
semiclassical energies or resonances are the poles of the 
response function g sc (E) in Eq. Q. 

The hydrogen atom in a magnetic held possesses a scal- 
ing property, as already introduced in Sec. [TTJ i.e., the 
classical dynamics does not depend separately on the en- 
ergy E and the magnetic held strength 7 but only on 
the scaled energy E = Ej~ 2 / 3 . Keeping the scaled en- 
ergy constant the response function can be written as a 
viz. 



function of w — 7 x ^ 



2 (w)=g S c{w)+Y,A vo e iws *° , 



(7) 
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where s po = 7 1 / 3 S l po is the scaled action of a periodic 
orbit and the amplitude 
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FIG. 8. (Color online) Number of periodic orbits with clas- 
sical action s po < s for all individual orbits (dashed green 
line), only one representative of a periodic-orbit bunch (dot- 
ted blue line), and one representative of a periodic-orbit bunch 
weighted with the size of the bunch (solid red line, nearly in- 
distinguishable from the dashed green line). All curves show 
an exponential increase, however, the number of periodic- 
orbit bunches grows with a significantly lower rate than the 
number of all orbits. 
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depends on the scaled action s ppo of the primitive peri- 
odic orbit (ppo) where the orbit is traversed only once, 
the repetition number r of that orbit, the leading eigen- 
value Ap po of the monodromy matrix, and the Maslov 
index ^ ppo . The Fourier transform of the scaled response 
function in Eq. ^ yields the semiclassical signal 



1 f°° 

C sc (s) = — J g sc (w)e~ lsw dw = ^„4 po r5(s - s po ) 



(9) 



which possesses <5-peaks with weight factors _4 po at the 
actions s — s po of the classical periodic orbits. 

The periodic orbit quantization is achieved by adjust- 
ing the semiclassical signal ([9| to its quantum analogue. 
The quantum resonances are the eigenvalues w = 7~ 1 / 3 
of the scaled Schrodinger equation at constant scaled en- 
ergy 



with resonances in the lower complex half plane, i.e., 
Im Wk < and residuals dk = 1 for non-degenerate states 
yields the quantum signal 
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(12) 



The semiclassical signal ([9| obtained from the classical 
periodic orbit data is now adapted to the quantum sig- 
nal (12 1, with the dk and Wk being free adjustable com- 



plex parameters by means of nonlinear signal processing 
methods. Technical details can be found in Refs. [22ll23| . 

The periodic orbit signal C sc (s) must be known in the 
interval < s < s max . The required signal length s max 
is proportional to the density of eigenstates Wk, which 
usually increases for higher excitations, i.e., an increasing 
signal length is required to resolve the denser parts of the 
quantum spectrum. 



B. Reduced data set 



-2E { f i 2 + v 2 ) + i(/iV+A 4 )-4 



= w~ 2 (A fi + A v )iJ with ^P=- p ^ p \P§p \ ■ ( 10 ) 

The Fourier transform of the quantum mechanical 
Green's function, 



ffqmH = 



w - w k 



(11) 



Since in classically chaotic systems the number of pe- 
riodic orbits grows exponentially with the orbital period 
or action, a linear increase of the length s max of the pe- 
riodic orbit signal ^ requires an exponential increase 
of the classical input data and the numerical effort for 
the periodic-orbit search. For the hydrogen atom in a 
magnetic field at scaled energy E = 0.5 the exponential 
increase of the number of trajectories with growing ac- 
tion is illustrated by the dashed green line in Fig. [8] For 
the construction of the periodic-orbit signal (J9J) with sig- 
nal length s ma x/27r = 17 more than one million periodic 



orbits must be computed. 

To reduce the required amount of periodic-orbit data 
we take advantage of properties of the periodic-orbit 
bunches. Since all orbits of a bunch are near action- 
degenerate (cf. Fig. [6]) and also have very similar stabil- 
ity properties and identical Maslov indices the parame- 
ters A po and s po in Eq. (|9| are approximately the same 
for all members of a periodic-orbit bunch. (Additional 
weight factors which arise from symmetry decomposition 
will be discussed in Sec. |TVC| ) The idea now is not to 
compute all members but only one representative (or very 
few representatives) of a periodic-orbit bunch, and to 
use the actions and appropriately weighted amplitudes of 
the representatives for the construction of the periodic- 
orbit signal. The Maslov index which determines the 
complex phase of the amplitudes in Eq. ([8|) can be ob- 
tained directly and without any numerical periodic orbit 
search from the symbolic dynamics of orbits and reads 
fipo = 3L — N + — N-, where L is the code length and 
N± are the numbers of plus and minus symbols in the 
code. Note that the Maslov index does not change by 
application of any of the four reconnection rules. 

The construction of the reduced data set is achieved 
as follows. We need the symbolic dynamics of all or- 
bits and the periodic-orbit parameters of representatives 
up to a certain action s max . In our computations we 
use Smax/27r = 20. In a first step we calculate the sym- 
bolic sequences of orbits. As some orbits with very long 
symbolic dynamics can contribute to the periodic-orbit 
signal the computing time is reduced by optimizing the 
order of generating the symbolic codes in a way that or- 
bits with short action are preferably obtained earlier than 
orbits with long action. The length limit of every sym- 
bolic sequence describing one periodic orbit is achieved 
if the approximate estimate of the action based upon the 
symbolic sequence of every orbit exceeds the action limit 
Smax- From the symbolic code of a periodic orbit its ac- 
tion can roughly be estimated by weighting the actions of 
the fundamental orbits 0, +, and - with the correspond- 
ing number of symbols in the code. We can calculate 
efficiently up to a maximum code length L = 16, i.e., the 
orbit data set is complete up to this length. Longer sym- 
bolic sequences are generated by adding minus symbols 
to shorter symbolic sequences. This does not change the 
shape of the orbits or the size of the groups in the most 
cases but adds only new loop parts in outward direction 
to the periodic orbits. In order to decrease computing 
time the same method is used to calculate all periodic 
orbits. 

We use the reconnection rules to group the periodic 
orbits. Arranging the orbits is a much faster procedure 
than the calculation of the actions because it is based 
only on integer arithmetic, which is very fast to handle on 
the computer. As the parameters of the periodic orbits 
in one group are very similar, we can choose one orbit as 
representative of the group and calculate all parameters 
of this trajectory. The results with appropriate weighting 
are then used for the construction of the periodic-orbit 



signal. 

The number of periodic-orbit bunches with action 
s po < s is shown as the dotted blue line in Fig. [8] Similar 
to the total number of orbits the number of the bunches 
also grows exponentially with the action, but with a sig- 
nificantly lower rate. For example at s/2n — 17 there are 
only about 80 000 bunches, compared to more than one 
million individual periodic orbits. 



C. Symmetry decomposition 

The diamagnetic hydrogen atom possesses a discrete 
symmetry viz. the reflection at the z — plane, and thus 
the quantum system has resonances in the decom- 
posed subspaces with even and odd z parity. In periodic- 
orbit theory the symmetry decomposition is achieved by 
multiplying the amplitudes A po in Eq. ^ with weight 
factors cr po , which depend on the chosen subspace and 
symmetry properties of the periodic orbits 24 j. The 
weight factors are 



(-1) 



for even parity, 
for odd parity, 



(13) 



with N+ the number of plus symbols in the symbolic 
code. In spectra with even parity all members of a 
periodic-orbit bunch have identical weight factors er po = 
1 and contribute with approximately the same ampli- 
tudes *Ap to the semiclassical signal in Eq. ffify, i.e., all 
orbits of a periodic-orbit bunch interfere constructively 
and the total weight of the bunch amplitude is approxi- 
mately the amplitude of a representative multiplied with 
the number of orbits in the bunch. In odd parity spectra 
the weight factors cr po of two orbits which are related by 
a single application of the reconnection rule R-l have dif- 
ferent sign and thus the contributions of the two orbits to 
the semiclassical signal in Eq. ^ approximately cancel. 
The total weight of the bunch is the nontrivial sum of the 
CT po of all individual orbits of the bunch, however, that 
weight can be obtained solely from the symbolic dynam- 
ics of orbits. 



D. Results and discussion 

For the hydrogen atom in a magnetic field exact quan- 
tum and semiclassical resonances have already been com- 
pared by Tanner et al. [20 . The quantum resonances 
were obtained by numerical diagonalization of the gener- 
alized eigenvalue problem ( 10 ) using a complex rotated 
complete basis set. The semiclassical resonances were 



computed with a modified and extended cycle expansion 
technique. As already mentioned at sufficiently high en- 
ergies E > E c = 0.329 the symbolic dynamics is com- 
plete without any pruning of orbits, however, trajectories 
going far away from the nucleus in the direction of the 
magnetic field axis were found to be marginally stable. 



9 





□ ' □ 


□ 




- a - 














• . 2S 




Q 


□ 






□ 


quantum mechanical values 




semiclassical reduced * 




semiclassical □ 





Re w 



(b) 



quantum mechanical values 
semiclassical reduced * 
semiclassical □ 



0.02 


-0.02 
-0.04 
-0.06 
-0.08 
-0.1 
-0.12 
-0.14 

I 

0.02 


-0.02 
-0.04 
I -0.06 
-0.08 
-0.1 
-0.12 

-0.14 I 1 1 

1 2 3 4 5 

Re w 

FIG. 9. (Color online) Resonances with (a) even and (b) 
odd parity of the diamagnetic hydrogen atom at scaled en- 
ergy E — 0.5. In most cases especially the real parts of the 
semiclassical resonances obtained with the reduced data set 
(red crosses) agree very well with the semiclassical resonances 
obtained with the complete periodic-orbit set (blue squares). 
For completeness and comparison the numerically exact quan- 
tum resonances are shown as black dots. 



Furthermore, resonances accumulate at the thresholds of 
the Landau channels. For these reasons the periodic or- 
bit quantization of the diamagnetic hydrogen atom is a 
very nontrivial and challenging task. 

Here we apply the harmonic inversion method intro- 
duced in Sec. |IV A| for period-orbit quantization. With 
a finite length signal this method can resolve sufficiently 
isolated resonances in the spectra but not the accumu- 
lation of resonances at the Landau thresholds. However, 
we do not focus on the comparison of the quantum and 
semiclassical resonances but on the results of the two 
semiclassical quantizations obtained with either the com- 
plete periodic-orbit set or with the reduced data set of 
the periodic-orbit bunches introduced in Sec. |IVB| 

Resonances with even and odd parity of the diamag- 
netic hydrogen atom at scaled energy E — 0.5 are pre- 
sented in Fig. [9j The semiclassical resonances have been 
obtained by harmonic inversion of the periodic-orbit sig- 



nal C sc (s) with s max /27r = 20 using either the complete 
periodic-orbit set (blue squares in Fig. [9]) or the reduced 
data set (red crosses). In the reduced data set we have 
used the two representatives with maximum or minimum 
action of each periodic-orbit bunch, which can be iden- 
tified by their symbolic codes. The numerically exact 
quantum resonances are shown by black dots in Fig. [9] 

As mentioned above the harmonic inversion of a finite 
length signal cannot reproduce resonances close to the 
accumulation points at, e.g., w = 1 and w = 3. For 
most of the semiclassical resonances with even parity in 
Fig. [9|a) the agreement is excellent for both the real and 
imaginary parts of the resonances. For the resonances 
with odd parity in Fig. |9^b) the agreement of the imag- 
inary parts is less perfect, however, the real parts still 
agree very well. The reason for the larger deviations 
in spectra with odd parity might be that the possible 
cancellation of periodic orbit contributions in this sub- 
space is more critical for the accuracy of the periodic- 
orbit signal than the solely constructive superposition of 
the contributions in the subspace with even parity. The 
results clearly demonstrate that the use of near action- 
degenerate periodic-orbit bunches allows for the reduc- 
tion of the classical data set by more than an order of 
magnitude and thus significantly increases the efficiency 
of periodic-orbit quantization methods. 



V. CONCLUSION 

The existence of periodic-orbit bunches in the classi- 
cally chaotic diamagnetic hydrogen atom has been re- 
vealed. The orbits of a bunch have a similar shape in the 
fundamental domain of the coordinate space and only 
differ in the behavior at self-encounters. We have intro- 
duced four reconnection rules which allow for the group- 
ing of orbits already on the level of the symbolic dynamics 
in terms of a ternary alphabet. 

The exponential proliferation of periodic orbits of 
a classically chaotic system implies that semiclassical 
methods based on Gutzwiller's trace formula typically 
require an exponentially growing classical data set to re- 
solve more eigenstates. We have shown that the use of 
one or few representatives of the near action-degenerate 
orbits can help to significantly improve the efficiency of 
semiclassical quantization methods. 

A peculiarity of the diamagnetic hydrogen atom is that 
periodic orbits at positive energies are marginally stable. 
It will thus be very interesting to investigate the exis- 
tence of near action-degenerate periodic orbits for other 
chaotic systems such as A^-disk or A-sphere scatterer, 
and to study the improvement of semiclassical quantiza- 
tion methods by using the periodic-orbit bunches. 



[1] M. C. Gutzwiller, Chaos in Classical and Quantum Me- [2] P. Cvitanovic and B. Eckhardt, 
chanics (Springer, New York, 1990). 



Phys. Rev. Lett. 



63, 823 



10 



(1989). 

J. Main, V. A. Mandelshtam, and H. S. Taylor, |Phys.| 
Rev. Lett.| 79, 825 (1997). 

J. Main, V. A. M andelshtam, G. Wunner, and H. S. Tay- 
lor, |NonnneSity[ll, 1015 (1998). 
J. Main, |Ph~ Rep. 316, 233 (1999). 
O. Bohigas, M. J. Giannoni, and C. Schmit, Phys. Rev 
52, 1 (1984). 



Lett. 



M. Sieber and K. Richter, Physica Scripta T90, 128 
(2001). 

M. Sieber, |jTPhys~A 35, L613 (2002). 
S. M Mer, S. Heusler, P. Braun, F. Haake, and A. Alt- 
land, |P^s7RevrTett] 93, 014103 (2004). 
S. Heusler, S. Miiller, P. Braun, and F. Haake, |J. Phys 
A 37, L31 (2004). 

S. M iiller, S. Heusl er, P. Braun, F. Haake, and A. Alt- 
land, [Phys~Revr]| 72, 046207 (2005). 
S. Heusle r, S. Miiller, A. Altland, P. Braun, and 
F. Haake, |Phys. Rev. Lettl 98, 044103 (2007). 
S. Miiller , S. Heusle~ A. Altland, P. Braun, and 
F. Haake, [New J. Phys.| ll, 103025 (2009). 
A. Altland, P. Braun, F. Haake, S. Heusler, G. Knieper, 



and S. Miiller, in Path Integrals - New Trends and 
Perspectives: Proc. 9th Int. Conf. (Dresden), edited by 
W. Janke and A. Pelster (World Scientific, Singapore 
2008) p. 40. 

[15] H. Hasegawa, M. Robnik, and G. Wunner, |Prog. Theor 

Phys SuppL] 98, 198 (1989). 

[16] H. Friedrich and D. Wintgen, Phys. Rep. 183, 37 (1989) 
[17] B. Eckhardt a nd D. Wintgen , | J. Phys~B| 23, 355 (1990) 
[18] K. T. Hansen, |Ph~ys. Rev. E| 51, 1838 (1995). 
[19] S. Schneider, diploma thesis, Ruhr-Universitat Bochum 

1997 (unpublished). 
[20] G. Tanner, K. T. Hansen, and J. Main, |Nonlinearity| 9 

1641 (1996). 

[21] J. Gehrke, diploma thesis, Universitat Stuttgart, 2008 
(unpublished). 

[22] J. Mai n, P. A. Dando, Dz. Belkic, and H. S. Taylor, J. 



Phys. A 33, 1247 (2000). 
[23] Dz. Belkic, P . A. Dando, J. Main, and H. S. Taylor, J. 
I Chem. Phys. 1 113, 6542 (2000). 

[24] P. Cvitanovic and B. Eckhardt, |Nonlinearity| 6, 277 
(1993). 



